Towards a Quantum Transport Description of Particles 
with finite Mass Width 1 



Jorn Knoll, GSI, Darmstadt 



J. KnollQgsi. de; http:/ /theory. gsi. dt 
Yu. B. Ivanov, Kurchatov Institute, Moscow 



D. N. Voskresensky, Moscow Institute for Physics and Engineering 

(Sept. 15, 1998) 



00 ■ 
ON 1 

On ; 

> : 
o . 

£: 

r- : 

<N . 
(N ' 

> : 

ON . 

■ 

On ; 
O . 
00 ■ 

5 : 
^ : 

a: 

Oh! 
D 

43 . 



X 



The effects of the propagation of particles which 
have a finite life time and an according width in their 
mass spectrum are discussed in the context of trans- 
port descriptions. In the first part the coupling of 
soft photon modes to a source of charged particles 
is studied in a classical model which can be solved 
completely in analytical terms. The solution corre- 
sponds to a re-summation of certain field theory di- 
agrams. The second part addresses the derivation of 
transport equations which also account for the damp- 
ing width of the particles. The ^-derivable method 
of Baym is used to derive a self-consistent and con- 
serving scheme. For this scheme a conserved energy- 
momentum tensor can be constructed. Furthermore, 
a kinetic entropy can be derived which besides the 
standard quasi-particle part also includes contribu- 
tions from fluctuation. 



I. INTRODUCTION 

With the aim to describe the collision of two 
nuclei at intermediate or even high energies one is 
confronted with the fact that the dynamics has 
to include particles like the delta or rho-meson 
resonances with life-times of less than 2 fm/c 
or equivalently with damping widths above 100 
MeV. The collision rates deduced from presently 
used transport codes are comparable in magni- 
tude, whereas typical mean kinetic energies as 
given by the temperature range between 70 to 
150 MeV depending on beam energy. Thus, the 
damping width of most of the constituents in the 
system can no longer be treated as a perturba- 
tion. 

As a consequence the mass spectrum of the 
particles in the dense matter is no longer a sharp 
delta function but rather acquires a width due to 
collisions and decays. The corresponding quan- 
tum propagators G (Green's functions) are no 
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longer the ones as in the standard text books for 
fixed mass, but rather have to be folded over a 
spectral function A(e,p), which takes a Lorentz 
shape A(e,p) ~ T/((e - e(P)) 2 + (r/2) 2 ) of width 
r/2 in simple approximations. One thus comes 
to a picture which unifies resonances which have 
already a width in vacuum due to decay modes 
with the "states" of particles in dense matter, 
which obtain a width due to collisions (colli- 
sional broadening). The theoretical concepts for 
a proper many body description in terms of a real 
time non equilibrium field theory have already 
been devised by Schwinger, Kadanoff, Baym and 
Keldysh [Q in the early sixties, extensions to rel- 
ativistic plasmas by Bezzerides and DuBois 
First investigations of the quantum effects on the 
Boltzmann collision term were given Danielewicz 
||, the principal conceptual problems on the 
level of quantum field theory were investigated 
by Landsmann 0] , while applications which seri- 
ously include the finite width of the particles in 
transport descriptions were carried out only in re- 
cent times, e.g. Pj5j-pd[. For resonances, e.g. the 
delta resonance, it was natural to consider broad 
mass distributions and ad hoc recipes have been 
invented to include this in transport simulation 
models. However, many of these recipes are not 
correct as they violate some basic principle like 
detailed balance j^j, and the description of reso- 
nances in dense matter has to be improved fiofl . 

In this talk the consequences of the propa- 
gation of particles with short life times is re- 
addressed and discussed. In the first part a 
genuine soft mode problem is studied: the cou- 
pling of a coherent classical field, the Maxwell 
field, to the stochastic Brownian motion of a 
charged particle. The rate of photons due to 
Bremsstrahlung, given by the classical current- 
current correlation function, can be obtained in 
closed analytical terms and discussed as a func- 
tion of the macroscopic transport properties, the 
friction and diffusion coefficient of the Brown- 
ian particle. The result corresponds to a partial 
re-summation of photon self energy diagrams in 
the real-time formulation of field theory. In the 
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second part of this talk a scheme is presented, 
how to come to a self-consistent, conserving and 
thermodynamically consistent transport descrip- 
tion of particles with finite mass width within 
the real-time formulation of non-equilibrium field 
theory. The derivation is based on and general- 
izes the ^-functional method of Baym . The 
first-order gradient approximation provides a set 
of coupled equations of time-irreversible gener- 
alized kinetic equations for the slowly varying 
space-time part of the phase-space distributions 
and retarded equations, which provides the fast 
micro-scale dynamics represented by the four- 
momentum part of the distributions. Functional 
methods permit to derive a conserved energy- 
momentum tensor which also includes corrections 
arising from fluctuations besides the standard 
quasi-particle terms. Memory effects [ fL3| appear- 
ing in collision term diagrams of higher order 
are discussed. The variational properties of <f>- 
functional permit to derive a generalized expres- 
sion for the non-equilibrium kinetic entropy flow, 
which includes corrections from fluctuations and 
memory effects. In special cases we demonstrate 
that the entropy can only increase with time (H- 
theorem). 

II. PRELIMINARIES 

The standard text-book transition rate in 
terms of Fermi's golden rule, e.g. for the photon 
radiation from some initial state \i) with occupa- 
tion rii to final states |/} 
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with occupation n w for the photon, is limited to 
the concept of asymptotic states. It is therefore 
inappropriate for problems which deal with par- 
ticles of finite life time. One rather has to go 
to the "closed" diagram picture, where the same 
rate emerges as 
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with now two types of vertices — and + for the 
time-ordered and the anti-time ordered part of 
the square of the amplitude. Together with the 
line sense and the — and + marks at the vertices 
a unique correspondence is provided between the 



oriented — > and — > propagator lines and the 
initial and final states. Thus such propagator 
lines define the densities of occupied states or 
those of available states, respectively. There- 
fore all standard diagrammatic rules can be used 
again. One simply has to extend those rules to 
the two types of vertices with marks — and + 
and the corresponding 4 propagators, the usual 

time-ordered propagator — > between two — ver- 
tices, the anti-time-ordered one ^+ between two 

+ vertices and the mixed — — ► or — 4 ones as den- 
sities of occupied and available states. For details 
I refer to the textbook of Lifshitz and Pitaevski 
(Til l . The advantage of the formulation in terms 
of "correlation" diagrams, which no longer refer 
to amplitudes but directly relate to physical ob- 
servables, like rates, is that now one is no longer 
restricted to the concept of asymptotic states. 
Rather all internal lines, also the ones which orig- 
inally correspond to the "in" or "out" states are 
now treated on equal footing. Therefore now one 
can deal with "states" which have a broad mass 
spectrum and which therefore appropriately ac- 
count for the damping of the particles. The cor- 
responding Wigner densities — — ► or — -» are then 
no longer on-shell 5-functions in energy (on-mass 
shell) but rather acquire a width in terms of the 
spectral function, e.g. for non-relativistic parti- 
cles 

iG-+ = ^ = Tf(p)A(p) 

iG+- = = (1 T f(p))A(p) (3) 

A{p) = T (P) 

(e + MF - A - Re E*(p)) * + (r(p)/2)2 

Here f(p) is the phase-space occupation at four- 
momentum p — (e,p), A is the spectral function 
with the damping width T(p) and in-medium on- 
shell energy 4 -Re E fl (p) and y, is the chemical 
potential. In general all quantities depend on 
both, energy e and momentum p. 

) » t 

t t+ co 

FIG. 1. Closed real-time contour with two exter- 
nal points x, y on the contour. 

The non-equilibrium theory can entirely be for- 
mulated on one special time contour, the so called 
closed time path fig. |l|, with the time argu- 
ment running from some initial time to to infin- 
ity and back with external points placed on this 
contour, e.g., for the four different components 
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of Green's functions or self energies. The special 
— h or H — components of the self energies define 
the gain and loss terms in transport problems, 
c.f. eq. |). 

III. BREMSSTRAHLUNG FROM 
CLASSICAL SOURCES 

For a clarification of the infra-red problem we 
first discuss two simple examples of soft modes 
in hard matter. These are examples in classi- 
cal electrodynamics, which both can be solved 
analytically to a certain extent: there the hard 
matter is described either by a diffusion process 
or by a random walk problem, respectively ]Tl)| . 
As the source particles move non-relativistically 
these cases do not suffer from standard patholo- 
gies encountered in the hard thermal loop (HTL) 
problem of QCD, namely the collinear singular- 
ities, where vq ks 1 and from diverging Bose- 
factors. The advantage of these examples is that 
damping can be fully included without violating 
current conserving and gauge invariants in the 
case of Abelian fields. The closed form results 
obtained correspond to a partial re-summation 
of certain planar diagrams, which just survive 
in the classical limit. The problem is related 
to the Landau-Pommeranchuk-Migdal effect of 
Bremsstrahlung in high energy scattering |lq ] . 
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FIG. 2. Current-current correlation function in 
units of e 2 < v 2 > as a function of time (in units 
of 1/T X ) for different values of the photon momen- 
tum q 2 = 3fc 2 r 2 /< v 2 > with k = 0, 1, 2, 3. 

The diffusion process is assumed to be de- 
scribed by a Fokker-Planck equation for the prob- 
ability distribution / of position x and velocity 
v 



Likewise fluctuations evolve in time by this equa- 
tion and this way determine the correlations. 
The two macroscopic parameters are the spatial 
diffusion coefficient D and a friction constant T x 
which determines the relaxation rates of veloci- 
ties (friction due to collisions with the medium). 
In the equilibrium limit (t — > oo) the distribution 
attains a Maxwell-Boltzmann velocity distribu- 
tion where T = m{V 2 )/3 = mDT x . The correla- 
tion function can be obtained in closed form and 
one can discuss the resulting time correlations of 
the current at different fixed values of the photon 
momentum q, fig. ^| (details are given in rcf. Jll|). 
For the transverse part of the correlation tensor 
this correlation decays exponentially as ~ e~ VscT 
at q — 0, and its width further decreases with in- 
creasing momentum q = Besides trivial kine- 
matical factors, the in-medium production rate 
is given by the time Fourier transform t — > uj. 
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FIG. 3. Rate of real photons d 2 N/(du>dt) in 
units of 4-7re 2 (« 2 ) /3 for a non-relativistic source for 
V x =50,100,150 MeV; for comparison the IQF results 
(dashed lines) are also shown. 

Fig. U displays the corresponding total pro- 
duction rates d 2 N / (dudt) of on-shell photons 
(number per time and energy; which is dimen- 
sionless) in units of 47re 2 (y 2 ) /3. One sees that 
the hard part of the spectrum behaves as ex- 
pected, namely, like in the IQF approximation 
the rate grows proportional to T x and this way 
proportional to the microscopic collision rate T 
(c.f. below). However independent of T x the 
rate saturates at a value of ~ 1/2 in these units 
around u> ~ T x , and the soft part shows the in- 
verse behavior. That is, with increasing colli- 
sion rate the production rate is more and more 
suppressed! This is in line with the picture that 
such photons cannot resolve the individual colli- 
sions any more. Since the soft part of the spec- 
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trum behaves like io/Y x , it shows a genuine non- 
perturbative feature which cannot be obtained 
by any power series in T x . For comparison: the 
dashed lines show the corresponding IQF yields, 
which agree with the correct rate for the hard 
part while they completely fail and diverge to- 
wards the soft end of the spectrum. For non- 
relativistic sources (t/ 2 ) <C 1 one can ignore the 
additional q-dependence (dipole approximation; 
c.f. fig. ||) and the entire spectrum is determined 
by one macroscopic scale, the relaxation rate T x . 
This scale provides a quenching factor 
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by which the IQF results have to be corrected in 
order to account for the finite collision time ef- 
fects in dense matter. 

In the microscopic Langevin picture one con- 
siders a classical process, where hard scatterings 
occur at random with a constant mean collision 
rate T. These scatterings consecutively change 
the velocity of a point charge from v m to v m +i 
to v m +2i ■ ■ ■ (in the following subscripts m and n 
refer to the collision sequence). In between scat- 
terings the charge moves freely. For such a mul- 
tiple collision process some explicit results can 
be given, since the correlated probability to find 
the charge at time t\ and t 2 at two different seg- 
ments with n scatterings in between follows from 
the iterative folding of the exponential decay law 
with decay time 1/r. Therefore the space inte- 
grated current-current correlation function takes 
a simple Poisson form 

m^oc J d 3 x 1 d 3 x 2 (j^(x 1 ,t - ^)r(x 2 ,t + ^)) 
= e 2 (v»(0)v v (t)) 
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with v = (l,v). This result represents a genuine 
multiple collision description of the correlation 
function. Here (•••)„> denotes the average over 
the discrete collision sequence {m}. This form, 
which one writes down intuitively, directly in- 
cludes what one calls damping in the correspond- 
ing quantum case. Fourier transformed it deter- 
mines the spectrum in completely regular terms 
(void of any infra-red singularities) where each 
term describes the interference of the photon be- 
ing emitted at a certain time or n collisions later. 

In special cases where velocity fluctuations are 
degraded by a constant fraction a in each col- 



one can re-sum the whole series in (g) and thus 
recover the relaxation result with 2T X (zr 2 ) = 
r ((v m — v m +i) 2 ) at least for q = and the cor- 
responding quenching factor (|5|). 

This clarifies that the diffusion result repre- 
sents a re-summation of the Langevin multiple 
collision picture and altogether only macroscopic 
scales are relevant for the form of the spec- 
trum and not the details of the microscopic col- 
lisions. Note also that the classical results, both 
for the diffusion equation (c.f. fig. ||) and for 
the Langevin process fulfill the classical version 
(h — > 0) of the sum rules discussed in refs. M,pT|. 



IV. RADIATION ON THE QUANTUM 
LEVEL 

We have seen that on the classical level the 
problem of radiation from dense matter can be 
solved quite naturally and completely at least 
for simple examples, and figs. || and |3| display 
the main physics. They show, that the damping 
of the particles due to scattering is an important 
feature, which in particular has to be included 
right from the onset. This does not only assure 
results which no longer diverge, but also provides 
a systematic and convergent scheme. On the 
quantum level such problems requires techniques 
beyond the standard repertoire of perturbation 
theory or the quasi-particle approximation. 




lision, such that (v ri 



The production or absorption 
rates are given by photon self + 
energy diagrams of the type 
to the right with an in- and 
outgoing photon line (dashed) . The hatched loop 
area denotes all strong interactions of the source. 
The latter give rise to a whole series of diagrams. 
As mentioned, for the particles of the source, e.g. 
the nucleons, one has to re-sum Dyson's equation 
with the corresponding full complex self energy 
in order to determine the full Green's functions 
in dense matter. Once one has these Green's 
functions together with the interaction vertices 
at hand one could in principle calculate the re- 
quired diagrams. However both, the computa- 
tional effort to calculate a single diagram and the 
number of diagrams, are increasing dramatically 
with the loop order of the diagrams, such that 
in practice only lowest order loop diagrams can 
be considered in the full quantum case. In cer- 
tain limits some diagrams drop out. We could 
show that in the classical limit of the quantum 
description, which in this case implies the hierar- 
chy oj, \q\, r <C T <C m together with low phase- 
space occupations for the source, i.e. f(x,p) <C 1, 
only the following set of diagrams survive 
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In these "Langevin" diagrams the bold lines de- 
note the full nucleon Green's functions which also 
include the damping width, the black blocks rep- 
resent the effective nucleon-nucleon interaction in 
matter, and the full dots the coupling vertex to 
the photon. Each of these diagrams with n inter- 
action loop insertions just corresponds to the n th 
term in the classical Langevin result (H). Thus 
the classical multiple collision example provides 
a quite intuitive picture about such diagrams. 
Thereby the diagram of order n describes the 
interference of the amplitude where the photon 
is "emitted" at some time and that where it is 
"emitted" n collisions later. Further details are 
given in pT|. 



V. ^-DERIVABLE APPROXIMATIONS 

Following Luttinger, Ward and the re- 

formulation by Cornwall, Jackiw and Tomboulis 
jl8| using path-integral methods for equilibrium 
case, the generating functional T{cj>, G} for the 
equations of motions, both, for the classical fields 
cf> =< <f> > and Dyson's equation for the propaga- 
tors G, can be expressed in terms of an auxiliary 
functional <&, where <!> is solely given in terms of 
full, i.e. re-summed, propagators G and full clas- 
sical fields 4>. Following p9[ it can be generalized 
to the real time case with the diagrammatic rep- 
resentation 

ir{0,G} = ir°{G°}+ <j> dxC ^^^} 
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E- 
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Here upper signs relate to fermion quanti- 
ties, whereas lower signs, to boson quantities. 
Thereby counts the number of self-energy 
^-insertions in the ring diagrams, while for the 
closed diagram of $ the value n\ counts the num- 
ber of vertices building up the functional <F Due 
to this factor such a set of diagrams is not resum- 
able in the standard diagrammatic sense. The r° 
solely depends on the unperturbed propagator 
G° (thin line) and, hence, is treated as a con- 
stant with respect to the functional variations 
in G(x,y) and 4>(x). The diagrams contribut- 
ing to <!> are given in terms of full propagators G 
(thick lines) and classical fields (p. As a conse- 
quence, these ^-diagrams have to be two-particle 
irreducible (label c2), i.e. they cannot be decom- 
posed into two pieces by cutting two propagator 
lines. The latter property matches diagrammatic 
rules for the re-summed self-energy H(x,y) and 
the source current J(x) of the classical field equa- 
tions, which results from functional variation of 
$ with respect to any propagator G(y,x), i.e. 



iS = =F<5i$/5iG, iJ = 6i<f>/5(/). 



(9) 



i$ {</>, G} 



It directly follows from the stationarity condition 
of r (||) with respect to variations of G and <fi on 
the contour 

ST {</>, G} /SG = 0, ST {4>, G} /5(j> = 0, (10) 

which indeed provides the Dyson equation with 
self-energy consistent with respect to the $- 
functional and the classical field equation. In 
graphical terms, the variation (^) with respect to 
G is realized by opening a propagator line in all 
diagrams of <&. The resulting set of thus opened 
diagrams must then be that of proper skeleton 
diagrams of S in terms of full propagators, i.e. 
void of any self-energy insertions. 

In order to arrive at a closed and consistent 
scheme we consider the so-called ^-derivable ap- 
proximation, first introduced by Baym | fl2"f based 
on ref. pp[ | within linear response to external per- 
turbation of equilibrated systems. They used the 
corresponding imaginary time formulation. A <f>- 
derivable approximation is constructed by con- 
fining the infinite set of diagrams for $ to ei- 
ther only a few of them or some sub-series of 
them. Note that $ itself is constructed in terms 
of "full" Green's functions and classical fields, 
where "full" now takes the sense of solving self- 
consistently the Dyson and Classical field equa- 
tion with the driving terms £ and J derived from 
this $ through relation (||). It means that even 
restricting ourselves to a single diagram in $, in 
fact, we deal with a whole sub-series of perturba- 
tion theory diagrams, and "full" takes the sense 
of the sum of this whole sub-series. Thus, a 
derivable approximation offers a natural way of 
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introducing closed, i.e. consistent approximation 
schemes based on summation of diagrammatic 
sub-series. In order to preserve the original sym- 
metry of the exact $ we postulate that the set of 
diagrams defining the ^-derivable approximation 
complies with all such symmetries. As a conse- 
quence, approximate forms of $( a PP r ) define ef- 
fective theories, where c&wpr.) serves as a g en _ 
erating functional for approximate self-energies 
j]( a PP r -) (x,y) and source currents J(x) through 
relation (|9|), which then enter as driving terms for 
the Dyson equations. The propagators solving 
this set of Dyson equations are still called "full" 
in the sense of the $( appr -)-derivable scheme. Be- 
low, we omit the superscript "appr." . 



VI. GENERALIZED KINETIC EQUATION 

A. Gradient Expansion Scheme 

For slightly inhomogeneous and slowly evolv- 
ing systems, the degrees of freedom can be subdi- 
vided into rapid and slow ones. Any kinetic ap- 
proximation is essentially based on this assump- 
tion. Then for any two-point function F(x,y), 
one separates the variable £ = (ii — £2,?'! — r3), 
which relates to rapid and short-ranged micro- 
scopic processes, and the variable X = i(t% + 
t2,fl + 7-2), which refers to slow and long-ranged 
collective motions. The Wigner transformation, 
i.e. the Fourier transformation in four-space dif- 
ference £ = x — y to four- momentum p of the con- 
tour decomposed components of F y ,i, j G { — h} 

F ij (X;p) = J dfe****' (X + £/2, X - f/2) 

(11) 

leads to a (co-variant) four phase-space formula- 
tion of two-point functions. The Wigner trans- 
formation of Dyson's equation ( |To| ) in { — +} no- 
tation is straight forward. For details and the ex- 
tensions to include the coupling to classical field 
equations we refer to ref. p9| . 

Standard transport descriptions usually in- 
volve two approximation steps: (i) the gradient 
expansion for the slow degrees of freedom, as well 
as (ii) the quasi-particle approximation for rapid 
ones. We intend to avoid the latter approxima- 
tion and will solely deal with the gradient approx- 
imation for slow collective motions by perform- 
ing the gradient expansion of the coupled Dyson 
equations. This step indeed preserves all the in- 
variances of the $ functional in a ^-derivable ap- 
proximation. 



B. Generalized Kinetic Equation in Physical 
Notation 

It is helpful to avoid all the imaginary factors 
inherent in the standard Green's function formu- 
lation and change to quantities which are real and 
in the homogeneous limit positive and therefore 
have a straight physical interpretation much like 
for the Boltzmann equation. We define 

F(X,p) = A(X,p)f(X,p) = i( T )G- + (X,p), 
F(X,p) = A(X,p)[l T f(X,p)} = iG + -{X,p) (12) 

for the generalized Wigner functions F and F 
and the corresponding four phase space distri- 
bution functions f(X,p) and Fermi/Bose factors 
[IT f(X, p)}. Here 

A{X,p) = — 21m G R (X,p) — F ± F (13) 

is the spectral function. According to retarded 
relations between Green's functions G 1 ^ , only two 
of these real functions are required for a complete 
description of the system 's evolution. 

The reduced gain and loss rates and total 
width of the collision integral are 

r in (x,p) = i( T )x- + (x,p), 
r aat (x, P ) = is+-{x, P ). (u) 

They determine the damping width 
r(A» = -2Im£ fl (A,p) 

= r out (A»±r in (A, P ), (15) 

where G R and T, R are the retarded propagator 
and self-energy, respectively. The opposite com- 
binations 

I(X,p) = T in (X,p) T r out (A,p), (16) 

determines the fluctuations. 

In terms of the new notation © - (0) and 111 
the first gradient approximation the generalized 
kinetic equation for F takes the form 

VF(X,p)-B = C(X,p) (17) 

with the differential drift operator (for simplicity 
in non-relativistic kinematics) 

„ / <9Re£ fl \ nu <9Re T, R d 
V ={ V »-^^) d x + ^X^8 p -» (18) 

with = (l,p/m). Further C(X,p) and 
B(X,p) are the collision and a fluctuation term, 
respectively 

C(X,p) = T ia (X,p)F(X,p) - T out (X,p)F(X,p) 
B = {T in ,RcG R }. (19) 



G 



We need still one more equation, which can 
be provided by the retarded Dyson equation. In 
terms of the new notation it takes the simple form 

VG R (X,p) + i {T,G R } =0, (20) 



M(X,p) + ^T(X,p))G R (X,p) = l, (21) 



with the "mass" function 

M(X,p)=p -^f-ReX R (X,p), (22) 

which relates to the drift operator via Df = 
{M, f| for any four phase-space function /. Sub- 
set © - @ is solved by fjj 



G = 



M(X,p)+iT(X,p)/2 



(23) 



A(X,p) 
Re G R (X,p) 



M 2 {X,p)+T 2 {X 1 p)/A : 

M(X,p) 
M 2 (X,p)+T 2 (X,p)/A' 



The canonical equal-time (anti) commutation re- 
lations for (fermionic) bosonic field operators im- 
ply sum-rules for the spectral function which for 
non-relaticistic kinematics read 



dpp 
2tt 



A(X,p) = 1. 



(24) 



We now provide a physical interpretation of 
various terms in the generalized kinetic equa- 
tion (|17| ). The drift term VF on the l.h.s. of 
eq. ( Jirj ) is the usual kinetic drift term includ- 
ing the corrections from the self-consistent field 
Re E fl into the convective transfer of real and 
also virtual particles. In the collision- less case 
C = B = 0, i.e. VF = (Vlasov equation), the 
quasi-linear first order differential operator V de- 
fines characteristic curves. They are the standard 
classical paths in the Vlasov case. Thereby the 
four-phase-space probability F(X,p) is conserved 
along these paths. The formulation in terms of a 
Poisson bracket in four dimensions implies a gen- 
eralized Liouville theorem. In the collisional case 
both, the collision term C and the fluctuation 
term B change the phase-space probabilities of 
the " generalized" particles during their propaga- 
tion along the " generalized" classical paths given 
by V. We use the term "generalized" in order to 
emphasize that particles are no longer bound to 
their mass-shell, M = 0, during propagation due 
to the collision term, i.e. due decay, creation or 
scattering processes. 

The r.h.s. of eq. ( |l7| ) specifies the collision 
term C in terms of gain and loss terms, which also 



can account for multi-particle processes. Since F 
includes a factor A, C further deviates from the 
standard Boltzmann-type form in as much that 
it is multiplied by the spectral function A, which 
accounts for the finite width of the particles. 
The additional Poisson-bracket term 



B = {T in ,ReG R } = 
MY 

+ (M 2 + r 2 /4) 2 



m 2 - r 2 /4 
(M 2 + r 2 /4) 2 

{r m ,r} 



vr. 



(25) 



is special. It contains genuine contributions from 
the finite mass width of the particles and de- 
scribes the response of the surrounding matter 
due to fluctuations. This can be seen from the 
conservation laws discussed below. In particu- 



lar the first term in (25) gives rise to a back- flow 
component of the surrounding matter. It restores 
the Nocthcr currents as the conserved ones from 
the intuitively expected sum of convective cur- 
rents arising from the convective VF terms in 
(17). The second term of ( p5| ) gives no contribu- 
tion in the quasi-particle limit of small damping 
width limit and represents a specific off mass- 
shell response, c.f. p2p3|]. 



C. Conservations of the Current and 
Energy— Momentum 

Special combinations of the transport equa- 
tions (|l7]) and the corresponding one for F 
weighted with e and p v , and integrated over 
momentum give rise to the charge and energy- 
momentum conservation laws, respectively, with 
the Noether charge current and Noether energy- 
momentum tensor defined by the following ex- 
pressions 



j»(X) = pr J (p(x,p) T F(X,p) 

&»{X) = \tt J (F(X,p) T F(X,p) 

(£ iat (X) - £ pot (X)) . (26) 



Here 



£' int (X) 



C im (X)) = 



5$ 



5\{x) 



(27) 



is the interaction energy density, which in terms 
of <I> is given by a functional variation with 
respect to a space-time dependent coupling 
strength of £ int -► X(x)Z int , c.f. ref. jT|. The 
potential energy density £ pot takes the form 
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where I = r in =p r out . Whereas the first 
term complies with quasi-particle expectations, 
namely mean potential times density, the sec- 
ond term displays the role of fluctuations / = 
Tin T Tout in the potential energy density. Since 
in many cases interaction and potential energy 
are proportional to each other, the same state- 
ment applies to the interaction energy, too. This 
fluctuation term precisely arises form the B-term 
in the kinetic eq. ( |l7|) , discussed around e g. ( |25| ) . 
It restores that the Noether expressions (|26| ) are 
indeed the conserved quantities. In this compen- 
sation we see the essential role of the fluctuation 
term in the generalized kinetic equation. Drop- 
ping or approximating this term would spoil the 
conservation laws. Indeed, both expressions in 
( ^6| ) comply exactly with the generalized kinetic 
equation (|17|), i.e. they are exact integrals of the 
generalized kinetic equations of motion. As usual 
the existence of such conservation laws require 
certain invariances which lead to certain consis- 
tency relations. In ref. jl9],^3| it has been shown 
that these are met if all the self-energies are $- 
derivable. 

In the field theoretical case there are contri- 
butions in (p6|), describing modifications of the 
vacuum-polarization in matter. These terms are 
generally ultra-violet divergent, and hence, have 
to be properly renormalized on the vacuum level. 
Alongside the spectral sum-rule ( gjj ) gets modi- 
fied by wave-function renormalization. 



D. Multiprocess Decomposition of 
$-Derivable Collision Term 

To be specific we consider a system of fermions 
interacting via a two-body potential V — Vq8{x — 
y), and, for the sake of simplicity, disregard 
its spin structure, by reducing spin and anti- 
symmetrization effects to a degeneracy factor d. 
To derive the decomposition of a ^-derivable col- 
lision term, we employ the same rules as de- 
scribed in ref. p3| . 

In the first example we consider the generating 
functional $ to be approximated by the following 
two diagrams 



i$ 




(28) 



the dashed line illustrating the decomposition. In 
the { — h} matrix notation of the Green's func- 
tions one can easily see that one-point diagrams 
do not contribute to the collision term, while de- 
composing the second one along the dashed line 
leads to a purely local result 



(2) 



xS 4 (p- 



d 4 pi d 4 p 2 d 4 p 3 



X 



(2tt) 4 (2tt) 4 (2tt) 4 
P1-P2- Pa) foFzFF! - F 2 F 3 FF 1 



(29) 



where the brief notation F t = F(X 1 p i ) etc. is 
used. This collision integral has precisely the 
form of the binary collision term of Boltzmann- 
Uchling-Uhlcnbcck (BUU), except for the fact 
that the distribution functions are not con- 
strained by the mass shell. The binary transition 
rate 



4 2) = vt = 



X 



(30) 



is non-negative in this case. 

The picture becomes more complicated, if $ 
involves diagrams of higher orders. For instance, 
let us add the following three point diagram to 
$, which is next in a series of ring diagrams, i.e. 



i$ = i ($(! 



(i) 



(2) 



(3) J 




(31) 



where one possible decomposition is illustrated 
by the dashed line. The corresponding self- 
energy becomes 




Now the collision term contains a non-local part 
due to the last diagram. The local part can easily 
be derived in the form 



C, 



loc 
(2) 



C(I) =d 2 



d 4 pi d 4 p 2 d 4 p 3 
(2tt) 4 (2tt) 4 (2tt) 4 



xS 4 (p + Pl - p 2 - p 3 ) {f 2 F 3 FF x - F 2 F 3 FF^ , 

where all the vertices in the off-shell scattering 
amplitudes are of the same sign, say " — " for 
definiteness, i.e., there are no " H — " and " — h" 
Green's functions left. The quantity + Cfe 
is again of the Boltzmann form 



R 



(2) 



Ri 3) = 



X 



IMS 
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the sub-label 2 denoting that 2 pairs of particle- 
hole lines are affected by the decomposition cut. 
It can be shown that under normal circumstances 
also this rate coefficient is positive. 



E. Kinetic Entropy 

Ignoring higher order gradients the generalized 
kinetic equation ( [l7|) provides us with the follow- 
ing relation 

V(aO = £ J ^Lln^C a (x,p), (34) 



where the quantity 




d 4 p 
(2^ 

TF a In - F a In 



ReG f ( T ^ln^-S^ln Fa 



A a dpu A a 



, (35) 



obtained from the l.h.s. of the kinetic equation is 
interpreted as the Markovian part of the entropy 
flow. Here we have introduced a summation over 
a denoting the different particle species and in- 
trinsic quantum numbers for a multi-component 
system. The interesting aspect is that for spe- 
cial local collision terms C a as the ones discussed 
above the r.h.s. of (Q) can be shown to be 
non-negative and therefore gives rise to an H- 
theorem. Again the functional properties of <f> 
have been used. The positivity of r.h.s. of ([34] is 
exactly given for $-functionals with two internal 
points for which in the equilibrium limit the zero 
component of the non-equilibrium entropy flow 
( |35| ) agrees with the corresponding equilibrium 
entropy. Memory corrections as contained in $- 
functionals with more than two points give rise 
to extra gradient terms which contribute to the 
entropy flow. For details we refer to our forth- 
coming paper p3|. 



VII. CONCLUSION 

In the first part of this talk the problem of soft 
modes in hard dense matter is discussed under 
circumstances which can be treated completely 
in analytical terms. The hard modes are de- 
scribed by a Fokker-Planck equation. They cou- 
ple to a classical Maxwell field for the soft modes, 
c.f. Fig. [|. For this Abelian case the result is 
conserving and completely gauge invariant even 



though the damping of the source particles is 
fully included. The friction coefficient (closely re- 
lated to the damping) determines the scale that 
separates soft from hard modes. This classical 
scheme is seen to re-sum a certain set of planar 
diagrams in the quantum case, which do survive 
in the classical limit. Such concepts are quite a 
general. In recent times they have been applied 
to the hard thermal loop (HTL) re-summation 
|2~4|] in terms of classical transport pslftq] . In 
the non-abelian QCD case, however, in order to 
preserve gauge invariance, the transport part is 
limited to the collision-less Vlasov equation, i.e. 
neglecting the damping of the source particles. 
A historical hard-loop re-summation scheme is 
the Fermi-liquid problem, where soft RPA-modes 
are treated by the coupling to the fermions in 
the Fermi-sea, the latter representing the hard 
modes, c.f. Fig. 0. 



QED: 



QCD: 



Fermi 
Liquids: 



Maxwell Fields 



Fokker Planck 
Eq. 



CI. Yang-Mills 
Fields 



Vlasov Eq. 
no Coll. -term 



RPA-modes 



Fermi-sea; -Bp 



FIG. 4. Hard Loop Re-Summation 

In practical terms we have seen that the spec- 
trum of soft particles resulting from collisions 
in dense matter can no longer appropriately be 
described by the quasi-particle approximation, 
since it leads to divergent results in the soft limit. 
Rather the finite time between successive colli- 
sions and the ensuing relaxation rates T x in dense 
matter lead to a considerable quenching of the 
production rate, e.g. at small photon energies. 
This can be compiled in the simple quenching 
factor (|^) . Fig. |3| summarizes the main behavior, 
which also is relevant in a quantum treatment of 
the source. 

In the second part a scheme is presented 
that leads to self-consistent conserving transport 
equations. There we essentially followed ideas 
suggested by Kadanoff and Baym in particular. 
The central quantity is a functional $ which gen- 
erates the driving terms for the classical field and 
transport equations. It can be truncated at any 
desired loop order and still provides equations 
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which fulfill conservation laws. We explicitly con- 
structed the energy momentum tensor for this 
derivable method. The gradient approximation 
provided equations of classical type for the phase- 
space distribution functions in four dimensions. 
At no place the quasi-particle approximation was 
necessary. Alongside from the <I>-derivable prop- 
erties a kinetic entropy could be derived, which 
in some cases leads to an H-theorem. 

In summary the method has the following ad- 
vantages: 

X provides a self-consistent & conserving 

transport scheme; 
X allows to include classical fields (soft 

modes); 

X includes all QM effects that are accounted 
for in the corresponding equilibrium treat- 
ment; 

X has no limitation to small widths; 
X includes delay-time, drag & back flow, and 
memory effects. 

There are two limitations: first, the derivation 
is limited to slow space-time variations of the 
macroscopic quantities; secondly, local symme- 
tries, like gauge invariance, may be violated by 
such re-summation schemes. The latter problem 
is inherent to all approaches, based on truncated 
self-consistent dynamical equations. 

Our considerations are of particular impor- 
tance for the theoretical description of nucleus- 
nucleus collisions at intermediate to relativistic 
energies. The kinematical feature are such that 
damping effects play an essential role, i.e. the 
energy uncertainty of the particles is compara- 
ble with their mean kinetic energy! In particular 
the bulk production and absorption rates of all 
particles with masses less than T, if calculated in 
standard quasi-particle approximation, are seri- 
ously subjected to the here discussed effects. 

In summary, the combined effort from many 
sides to include the finite width of the parti- 
cles in dense matter, may give hope for a unified 
transport theory which appropriately describes 
both, the propagation of resonances and of off- 
shell particles in the dense matter environment. 
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